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ABSTRACT 


This  report  contains  a comprehensive  survey  of  theoretical 
work  on  the  fluid  dynamic  effects  of  roughness  in  lubrication. 

Of  necessity,  much  of  the  report  deals  with  Reynolds 
roughness;  that  is,  roughness  which  can  be  locally  analyzed 
with  Reynolds  equation  because  the  wavelength  of  the  roughness 
is  not  too  short.  One-sided  and  two-sided  striated  roughnesses 
are  considered.  Statistical  and  multiple-scale  treatments  are 
discussed  and  intercompared.  Finally,  the  effects  of  short- 
wavelength  roughness  (Stokes  roughness)  are  treated,  as  well 
as  the  consequences  of  slip  in  highly  rarefied  situations. 

Some  new  results  for  isotropic  Reynolds  roughness,  for 
Stokes  roughness,  and  for  molecular  slip  are  given  preliminary 
exposure. 
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A REVIEVJ  OF  THEORIES  FOR  THE 


FLUID  DYNAMIC  EFFECTS  OF  ROUGHNESS 
ON  LAMINAR  LUBRICATING  FILMS 


H.G.  Elrod 


I.  Introduction; 

The  purpose  of  this  paper  is  to  summarize  some  theoret- 
ical results  concerning  the  effects  of  surface  roughness  on 
laminar  lubricating  films.  Old  and  new  work  are  inter-relat- 
ed, in  an  attempt  to  show  what  is  known,  and  what  needs  to  be 
found  out.  For  potential  investigators,  a considerable  bib- 
liography is  appended. 

Surface  roughness  may  be  deliberately  created,  as  in 
grooved  bearings,  to  promote  a pumping  effect.  It  may  be  the 
consequence  of  a fabrication  process,  such  as  drilling,  broach- 
ing, grinding,  lapping.  Or  it  may  be  the  consequence  of 
scratches  produced  by  debris.  Obviously,  then,  striations  will 
be  associated  with  many  forms  of  surface  finish. 

For  present  purposes,  it  suffices  to  observe  that, 
quite  typically,  roughness  amplitude  excursions  have  associat- 
ed elevation  angles  of  10°.  Dependent  on  operating  conditions, 
this  amplitude  may  become  a substantial  fraction  of  the  film 
thickness, h, with  characteristics  wavelengths  ranging  from 
0 (h)  -•  «. 

In  terms  of  film-thickness  reduction,  the  analyses  to 
be  discussed  herein  are  to  be  regarded  as  the  last  purely 
fluid  dynamical  treatments  before  some  solid-to-solid  inter- 
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action  of  asperities  takes  place.  The  flow  is  treated  as  lam- 
inar and  Newtonian,  with  either  liquid  or  gas,  unless  other- 
wise specified.  Cavitation  or  film  rupture  sometimes  occurs 
in  conjunction  with  roughness,  as  in  face  seals.  But  these 
effects  are  not  examined  here.  On  the  other  hand,  molecular 
slip,  which  takes  place  in  ultra-thin  gas  passages,  is  dis- 
cussed . 


Figure  1.1  Striated  Roughness 

Figure  1.1  shows  the  striated  roughness  with  which 
this  paper  is  mostly  concerned.  The  film  thickness  is  a 
function  only  of  the  coordinate,  z.  Its  cross-section  may 
be  arbitrary,  but  is  almost  periodic.  The  moving  surface  is 
smooth,  and  has,  by  choice  of  angle,  "B",  only  an  x-wise  com- 
ponent, U.  When  B=tt/2,  the  striations  are  called  "transverse", 
and  when  B=0,  the  striations  are  "longitudinal". 

If  the  characteristic  wavelength,  A,  of  the  striations 
is  large  compared  with  the  gap,  h,  then  it  is  legitimate  to 
use  Reynolds  equation  to  evaluate  roughness  effects.  Such 
roughness  we  term  "Reynolds  Roughness".  On  the  other  hand, 
if  the  wavelength  is  so  small  as  to  necessitate  the  use  of  the 
full  Stokes  equation,  we  use  the  name  "Stokes  Roughness". 
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II.  Analysis  of  Discrete  Rippling. 

The  earliest  analyses  of  roughness  effects  were  made 
by  means  of  Reynolds  equation  for  simple  cross-sections  si- 
nusoidal, rectangular,  triangular) . Michell  (1)  and  Burton 
(4)  used  sinusoidal,  first-order  linearization  treatments, 
the  former  with  longitudinal  striations  on  a simple  slider 
bearing,  and  the  latter  with  transverse  striations  in  a par- 
allel-plate thrust  situation.  Such  analyses  are  useful  in 
•disclosing  some  possible  effects  of  roughness,  but  do  not 
yield  the  load,  which  is  a second-order  effect.* 

Finite-difference  solutions  of  Reynolds  equation  in 
the  presence  of  roughness  were  obtained  by  Cameron  (3) (trans- 
verse) , Dowson  & Whomes  (7) (longitudinal) , Shelly  and  Ettles 

(8)  (longitudinal  & transverse) , and  by  Tender  & Christensen 

(9)  (longitudinal) . Because  of  the  necessity  of  numerical 
description  of  individual  grooves,  it  is  not  feasible  to  in- 
corporate in  analyses  of  finite  bearings  a large  number  of 
striations.  Thus  the  hypothesized  surfaces,  although  retain- 
ing validity  in  their  own  right,  are  not  realistic  renditions 
of  typical  roughness  profiles  ....  their  roughness  wavelength 
is  too  long.  Nevertheless,  some  of  these  numerical  studies 
show  the  important  fact  that  the  characteristic  wavelength  of 
the  roughness  ceases  to  affect  Reynolds-equation  solutions, 
even  when  it  is  only  moderately  small  compared  with  overall 
bearing  dimensions. 


* Footnote:  Although  reference  (1)  is  widely  quoted  as 
providing  the  effects  on  load  due  to  longitudinal  roughness, 
the  small-roughness  portion  of  the  analysis  actually  just  ex- 
hibits the  effect  of  enlarged  mean  clearance. 
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BIA  = 0 37. 

• Points  of  inflection. 


Figure  2.1  Assumed  Waveform  of 
Machining  Marks 
(from  ref.  2.8) 

Figure  2.1  shows  the  profile  assumed  by  Shelly  and 
Ettles  to  mimic  machining  marks  from  aero  boring  (longitudi- 
nal) or  from  broaching  (transverse) . They  studied  the  effects 
on  the  performance  of  a 50.8  mm  diam.  journal  bearing,  17.8  mm 
long,  with  incompressible  lubricant. 

Figure  2.2  presents  their  results  for  a protuberance 
amplitude  (h^)  of  1.6  ym  in  a clearance  (c)  of  25.4  pm.  The 
roughness  wavelength  was  varied  to  produce  different  numbers, 
n,  of  waves  either  circumferentially  (longitudinal)  or  axially 
(transverse) . Solutions  include  cavitation  with  Reynolds 
boundary-condition  . Note  the  dramatically  rapid 

approach  of  the  load  capacity  to  an  asymptotic  value.  (The 
indicated  load  reductions  are  the  consequence  of  basing  the 
reference  "smooth”  loads  on  minimum,  rather  than  mean,  film 
thickness) . The  early  work  of  Cameron  (3)  with  transverse 
waves  on  an  inf initely-wide  slider  (see  Fig. 2. 2),  and  the  more 
recent  work  of  Tender  and  Christensen  (9)  with  longitudial 
waves  on  cylindrical  rollers  exhibit  the  same  asymptotic  effect. 


Direction  Forward  ThruM 

of  wave*  motion  bearing 

Longinidintl  Rotation  Infinity  long 

iquecze  (amplitudc/A.„in) 
- 0 67 

Traniverac  Roution 
tqueeze 

Cameron 


Figure  2.2 

Asymptotic  Effect 
of  Number  of  Waves, 
(from  ref. 2. 8) 


For  the  special  case  of  roughness  in  the  form  of  rec- 

tilinearly  striated  rectangular  grooving,  Reynolds  equation 

2 

reduces  to  Laplace's  eq.,  *7?  = 0, 

in  each  of  the  constant  film-thicknesses  regions.  Muijderman 
(10)  showed  an  exact  solution  could  be  obtained  for  the  pres- 
sure distribution  by  the  juncture  of  linear  patterns,  (see 
Fig.  2.3)  all  conditions  being  satisfied  except  for  edge  ef- 
fects at  the  entrance  & exit  of  the  grooving.  In  the  case 
shown,  there  are  pressure  ripples  across  the  grooving  super- 
imposed an  a general  buildup  along  the  grooving.  Subject  to 
the  limitations  of  Reynolds  equation  and  the  hypothesis  of  in- 
compressible Newtonian  lubricant,  the  excursions  of  the  pres- 
sure ripple  can  be  very  substantial/  and  need  only  be  small 
compared  with  the  total  pressure  buildup  from  the  pumping  ac- 
tion/ to  minimize  edge  corrections. 


Figure  2.3  Pressure  Ripples  Generated 

Across  Rectangular  Roughness 
Elements  (from  ref. 10) 
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III,  Averaged  Variables  with  One-Sided  Reynolds  Roughness.* 

Even  with  the  deliberate  roughness  of  herringbone  and 
spiral-groove  bearings,  where  a/L  is  much  greater  than  for 
typical  surface  roughness,  it  is  desirable,  if  possible,  to 
deal  with  average  pressures  and  average  film  thicknesses.  On 
good  aul:hority,  the  writer  has  been  told  that  an  analysis 
achieving  this  purpose  was  accomplished  by  Boeker  circa  1944 
(11a).  However,  the  first  publicly  disseminated  analysis  is 
due  to  Whipple  (31b). 

The  analysis  of  Whipple, as  well  as  those  of  Muijderman 
(10),  Hirs  (12),  and  Elrod  (13),  takes  advantage  of  the  pres- 
sure periodicity  associated  with  parallel  opposing  surfaces, 
with  surfaces  uniformly  grooved.  Whipple  and  Muijderman  con- 
sider rectangular  grooving  only,  whereas  Hirs  and  Elrod  ac- 
count for  "arbitrary"  roughness.  We  suimnarize  briefly  here 
Hirs'  analysis. 

Figure  1.1  illustrates  the  situation  analyzed.  The 
smooth  upper  face  (towards  reader)  moves  with  velocity  U in 
the  x-direction.  Pressures  at  y=0  and  y=L  are  uniform  at  Pq 
and  Pj^,  respectively.  The  bearing  has  infinite  extent  in  the 
x-direction.  Without  proof,  Hirs  makes  the  following  plausi- 
ble statements  for  the  narrow-groove  limit  (i.e.,  when  A/L<<1) 
with  an  incompressible  fluid. 

a)  4^  is  independent  of  position  (an  hypothesis  made 

0 S 

earlier  by  Michell(2.1) 

b)  , the  film  flow  in  the  z-direction,  is  independent 
of  position. 


♦ The  roughness  pattern  must  not  move  with  respect  to 
the  gross  mean-film  configuration. 
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Applying  Reynolds  relations  for  the  volumetric  flows, 
one  obtains  from  (a) : 


12  m 3s  m 


■^■U  (cosB)  (h) 

2 m 


[3.1] 


where  the  subscript  "m"  indicates  an  average  over  one  wave- 
length, A,  that  includes  one  groove-ridge  pair. 


In  addition,  one  obtains  from  (b) : 


— } " ■ ~ ■jU(sinB){-^^ — } 

^ "*Mh3)-'m  ™ ^ Mh2)J 


m 


[3.2] 


Since  spatial  differentiation  and  integration  commute: 


9£.  ^ 3£ 


3s  m 


(■^) 
'3s  m 


(12)  = (12) 

3z  m 3z  m 


[3.3] 


where  p is  the  local  mean  pressure.  + With  and  con- 
stant, both  eqs . 3.1  & 3.2  can  be  integrated  subject  to  pres- 
sure or  flow  boundary  conditions.  The  performance  of  a simple 
herringbone  bearing  can  thereby  by  computed. 

Details  have  been  given  here  because  the  same  equations 
are  crucial  to  the  theory  of  Christensen  (see  section  IV) , and 
because  it  is  feasible,  even  in  a review  paper,  to  do  so.  We 
shall  reserve  comment  on  the  assumptions  (a)  and  (b)  until 
section  V. 

In  ref,  (13)  the  same  problem  is  attacked  by  assuming 
that  the  pressure  can  be  expanded  in  terms  of  n^y/L  and  i;=z/A 
with  t:A/L  as  a small  parameter.  Thus: 


+ Strictly,  the  p's  for  the  s-  and  z derivatives  arc 
not  the  same  except  for  linear  pressure  profiles,  but  the 
discrepancy  vanishes  as  A/L-»0 . 
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P = PQ(n,c)  + Ep^(r),c)  + E2p2(n,c)  etc.  [3.4] 

This  expansion  is  used  in  connection  with 

h = hQ  + Ef(0  [3.5] 

and  Reynolds  equation,  and  periodicity  is  invoked.  Differen- 
tial equations  are  found  for  p^,  the  mean  pressure,  p^,  the 
first  order  pressure  ripple,  etc.  The  results  apply  for  in- 
compressible and  compressible  lubricants,  and,  in  a sense, 
legitimize  earlier,  more  intuitive  work. 

The  limitation  of  constant  mean  film  thickness,  embod- 
ied in  Whipple's  work,  was  removed  by  Vohr  and  Pan  (14),  again 
on  an  intuitive  basis. Their  formulations  are  for  the  rectangu- 
lar-grooving equivalents  of  eqs.  (3.1)  and  (3.2).  These  ex- 
pressions they  employ  with  local  mean  thicknesses  and  local 
mean  fluid  density,  p,  and  write: 

V.p3  = 0 [3.6] 

Many  useful  studies  of  grooved-bearing  performance  have  been 
made  employing  the  Vohr-Pan  differential  equation. 

In  ref.  (15)  a small-parameter,  multiple-scale  analy- 
sis is  used  to  generate  the  appropriate  differential  equations 
for  varying  mean  film  thickness, arbitrary  orientation  of  the 
roughness  striations,  and  arbitrary  roughness  cross-section. 
The  expansion  adopted  is: 

p = pQ(n,z,c)  + epj^(n,z,c)  + e2p2(ti,z,c)  etc.  [3.7] 

explicit  recognition  being  given  to  the  possible  of  slow  pres- 
sure variation  in  terms  of  the  "slow"  variable,  z,  and  fast 
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variation  (rippling)  in  terms  of  the  "fast"  variable,  c.  An 
analogy  might  be  the  expression  of  atmospheric  temperature  in 
in  terms  of  month  and  hour  of  the  day.  The  slow  seasonal 
temperature  trend  would  be  known  as  a function  of  the  month, 
and  the  daily  fluctuation  as  a function  of  the  hour.  Again, 
the  analysis  is  in  agreement  with  more  heuristic  formula- 
tions; namely,  eq.(3.6)  with  components  of  Q given  by  eqs. 
(3.1)  and  (3.3)  is  supported. 

Within  their  stated  regions  of  validity,  all  of  the 
foregoing  (section  III)  analyses  yield  results  identical  to 
those  to  be  given  in  the  next  section. 


IV.  Statistical  Treatments  of  One-Sided  Reynolds  Roughness. 


The  haphazard,  semi-random  distribution  of  roughness 
on  many  surfaces  has  lead  investigators  to  consider  statisti- 
cal methods  of  analysis.  The  first  such  treatment  was  made 
by  Tzeng  & Saibel  (16)  . Figure  (4.1)  shows  the  system 


Figure  4.1  Simple  Slider  Bearing  with 
Roughness  on  One  Surface 

analyzed.*  The  authors  write,  for  an  incompressible  fluid: 


dx 


[4.1J 


where  "h"  is  the  film  thickness  including  roughness.  One 
integration  yields: 


* Fig.  1 of  ref. (16)  shows  two-sided  roughness,  for 
which  their  analysis  is  not  applicable  because  of  the  neglect 
of  time-derivatives.  (See  section  VI)  The  same  comment  ap- 
plies to  refs.  (18,19,20). 


Now  the  authors  write: 


h = h(x)+0  [4.3] 

where  h is  the  mean  of  a population  of  bearings  with  similar 
nominal  specification  and  "e"  is  some  random  fluctuation  about 
the  mean,  being  different  from  bearing  to  bearing. 

Thus;  <h> (expected  value)  = <h+e>  = h [4.4] 

The  authors  now  hypothesize  that  the  through-flow,  Q, 
will  be  similar  for  all  bearings  of  this  type,  and  take  the 
expected  values  of  the  quantities  in  eq.(4.2).  The  result 
is : 


< 


dx  dx 


12(iQ<— > 

h3 


[4.5] 


The  boundary  conditions  on  <p>  are: 
<p(0)>  = <p(L)>  = 0 

Integration  of  eq.(4,5)  yields  the  load. 


Thus : 

L 

Load  = |<p>dx  [4.6] 

0 

To  use  these  equations,  some  knowledge  of  the  rough- 
ness distribution  is  necessary.  Tzeng  & Saibel  proposed  a 
so-called  "beta  distribution"  function  which  yields  finite 
maximum  deviations  from  a mean  asperity  height,  but  other- 
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wise  closely  resembles  a Gaussian  normal  distribution.  Chri- 
stensen (18)  followed  their  example  with  a slightly  different 
expression,  which  to  avoid  confusion,  is  the  only  one  to  be 
given  here: 


[4,7] 


According  to  this  formula,  the  fraction  of  asperities  at  a 
given  "x"  in  the  bearing  population  with  heights  (measured 
from  a mean)  between  h^^^  and  hg2  is: 

fSl 

f(e)de 


The  function  is  compared  with  the  normal  distribution  curve 
in  Fig. (4.2) 


Goufston  distributKy) 


Polynomiol  distribution 


.1  J I L 1— — - 

-5  -I  0 I 2 5 


Figure  4.2  Roughness  Distribution 
Curves  Compared 
(from  ref. 18) 


Note  that  the  maximum  amplitude  is  three  standard  deviations 
(i.e. , c=3a) . 
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Tzeng  & Saibel  cite  only  one  numerical  example,  but 
their  analysis  yields  the  curves  for  transverse  roughness 
shown  in  Fig. 4. 3. 


Figure  4.3  Effect  of  Roughness  on 

Load  Capacity  of  Infinitely 
Wide  Slider  Bearing 
(from  ref. 18) 

Christensen (18)  and  Christensen  & Tender (9 , 19 , 20)  adopt 
a somewhat  different  statistical  viewpoint.  They  regard  the 
function  f(e)  as  providing  the  roughness  distribution  for  the 
population  of  surface  elements  neighboring  a given  "x”  in  a 
given  bearing.  This  difference  in  approach  has  lead  to  mis- 
understanding , even  among  the  analysts  themselves (18) . 

Essentially,  the  Christensen  (18)  roughness  theory 
rests  on  the  same  two  postulates  adopted  by  Hirs  (12) , but 
with  the  Vohr-Pan  (14)  viewpoint  that  the  flow  relations  (3.1) 
and  (3.2)  are  locally  applicable  in  a bearing  with  variable 
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mean  film  thickness.  In  addition,  accenting  the  uncertainty 
of  the  actual  roughness  distribution,  they  interpret  all 
means  as  statistical  expectations. 

For  transverse  roughness,  6=0 (x,t)  or  6=^/2, eqs. (3.1) , 
(3.2) , and 

V • Q = — [4.8] 


leads  to  one  of  Christensen's  differential  equations: 


+ = ^ + 12y# 

3x  — =5  3x  ay  ay  ^ ax  — » ^at 

h"  h" 


[4.9] 


For  longitudinal  roughness,  0=6  (y,t)  or  6=0,  his  re- 
sult is:  + 

^ li  - ^ 4=  - 6vu||  ^ 12v||  14.101 

h"'’ 

Here  h^  stands  for  the  local  mean  value  or  expectation  of  h"^. 

In  addition,  for  uniform  isotropic  roughness  with 
steady-state  operation,  Christensen  & Tender  propose  (19,20): 


_ih7i£  + 
3x  ax 


» h?  |£  = 6„u|i 

ay  ay  9x 


[4.11] 


These  differential  equations  resemble  very  closely  the 
usual  Reynolds  equation,  although,  as  Christensen  points, 
there  is  no  completely  equivalent  smooth  bearing  to  be  found 
by  adopting  some  appropriately  chosen  mean  film  thickness. 

But  the  niamerical  treatments  are  almost  identical  and 

there  is  no  necessity,  either,  of  treating  every  little  rip- 
ple individually. 


* The  more  general  case  of  arbitrary  striation  angle, 
6,  and  variable  fluid  density  is  treated  in  (15). 
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Christensen  and  Tender  have  applied  differential 
equations  (4.9)  and  (4.10)  to  slider  and  journal  bearings  of 
finite  and  infinite  width  (18,19,20)  using  roughness  distri- 
bution function  (4.7).  Fig. (4.3)  shows  the  theoretical  load 
capacities  of  the  slider  bearing  in  Fig.  (4.1)  for  the  two 
types  of  roughness,  longitudinal  and  transverse.  It  is  seen 
that  the  longitudinal  roughness  has  scarcely  any  influence 
on  load,  whereas  the  transverse  roughness  certainly  can  pro- 
duce a measurable  increase  in  load  capacity  over  that  for  a 
perfectly  smooth  bearing.  Notice  that  the  curves  terminate 
where  c/mL=h/mL  since  there  is  then  incipient  contact  of  the 
asperities  with  the  smooth  opposing  surface. 

The  effect  of  transverse  roughness  on  slider  friction 
is  shown  in  Figures  (4.4)  and  (4.5). 


I 

I 

1 

Figure  4.4  Longitudinal  and  Transverse  ! 

Roughness  Effects  on  Slider 
Friction,  (from  ref. 16). 


From  Fig.  (4.4)  it  is  observed  that  longitudinal  roughness 
scarcely  affects  the  friction  at  all,  whereas  transverse 
roughness  can  increase  it  appreciably.  However,  for  a given 


This  is  perhaps  an  appropriate  place  to  summarize  ! 

some  of  the  different  viewpoints  embodied  in  the  analyses  of  | 

Sections  II, III  and  IV,  and  to  make  intercomparisons.  'i 

i ! 

ii 

In  Section  II  the  common  feature  is  the  assumption  of 
some  specific  roughness  contour  (rectangular,  sinusoidal, 
etc.)  to  which  Reynolds  equation  is  applied,  either  by  a | 

small-perturbation  analysis  or  by  a direct  numerical  method. 
Moreover,  the  detailed  interaction  between  the  surface  rough- 
ness and  the  lubricating  film  is  followed. 

In  Section  III,  on  grounds  of  periodicity  or  short 
roughness  wavelength,  detailed  interactions  are  "smeared" 
or  averaged,  and  differential  equations  are  developed  which  i 

permit  calcualtions  of  load  capacity,  friction  and  other  per- 
tinent bearing  quantities  for  wide  classes  of  bearing  config- 
urations. 

! 

In  Section  IV,  two  different  statistical  approaches  j 

are  used  to  derive  differential  equations  for  bearing  perfor- 
mance. The  second  of  these  approaches  has  now  been  quite  ex- 
tensively applied. 

I 

t ; 

We  have  already  seen  that  there  is  a high  degree  of 
overlap  in  applicability  between  the  most  general  of  the  de- 
velopments described  in  Section  II  (ref. 15)  and  the  most  gen- 
eral development  in  Section  III  (ref .18) (with  the  latter  hav- 
ing publication  priority) . We  have  also  seen  that,  for  rough- 
ness wavelengths  which  are  small  compared  with  overall  bearing 
dimensions,  the  results  of  Section  I are  in  accord  with  those 
of  II  and  III,  when  the  particular  form  of  roughness  is  taken 
into  consideration.  The  question  may  now  be  asked:  "To  what 
extent  does  the  particular  form  of  roughness  matter"? 


^ ^ 
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For  small-amplitude  roughness,  the  answer  is  given  in 
refs.  (13)  and  (18),  where  it  is  shown  that, in  this  case, the 
special  character  of  the  roughness  does  not  matter,  only  its 
rms  value.  Consider  the  differential  equations  (4.9)  and 
(4.10).  Four  means  of  the  film  thickness  appear  namely, 

h^  , h ^ , h ^/h  ^ and  h 

The  values  for  these  quantities  which  would  arise  for 
the  same  specified  h and  variance,  o,  are  now  in  question. 

If  all  roughness  shapes  resulted  in  exactly  the  same 
values  for  the  first  three  of  these  quantities,  identical 
pressure  distributions  would  be  predicted. 

Now  the  first  quantity,  h^  = (h+e)^  = h^+3a^ 
is  the  same  for  any  symmetrical  roughness  distribution.  Re- 

_3  _2  -3 

suits  for  the  other  quantities,  h and  h /h  are  compared 
in  Figs.  (5.1)  and  (5.2)  for  simple  sinusoidal  and  rectangu- 
lar roughness,  and  for  the  roughness  distribution  adopted  by 
Tzeng  & Saibel  and  Christensen  & T0nder  eq.(4.1).  Even  when 
o/h=0.2  (0.33  corresponds  to  incipient  solid-solid  contact), 

-3  -3  _2  - -3 

the  discrepancies  for  h h and  for  h /h  h do  not  exceed 
4%.  Therefore,  in  so  far  as  Reynold  roughness  analyses  are 
concerned,  there  is  virtually  nothing  to  choose  among  the 
different  forms. 

In  Section  IV,  the  difference  between  the  Tzeng  & Sai- 
bel and  Christensen  & Tender  statistical  approaches  was  de- 
fined. Why,  then,  do  the  authors  use  virtually  the  same  as- 
perity distribution  function?  Because,  as  a practical  mat- 
ter the  bearing  population  concept  (many  "like"  bearings)  is 
difficult  to  implement,  not  only  because  of  the  number  of 
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bearings  required,  but  because  of  the  virtual  impossibility 
of  ascertaining  positions  "x"  within  the  film  to  accuracies 
of  the  order  of  roughness  wavelength.  These  considerations 
do  not,  of  course,  invalidate  the  Tzeng  & Saibel  approach, nor 
their  slider  analysis.  In  actuality,  though,  Tzeng  & Saibel 
must  invoke  an  "ergodic  hypothesis",  namely  that  the  spatial 
distribution  of  asperities  on  a given  bearing  is  similar  to 
that  obtained  from  replication  of  the  bearing  sample.  This 
reconciliation  of  roughness  distributions  does  not,  however, 
render  the  difference  in  the  two  statistical  approaches  triv- 
ial, since  the  local  sampling  approach  of  Christensen  & Ten- 
der permits  them  to  generate  second-order  partial  differen- 
tial equations  applicable  to  a wide  variety  of  bearings. 

In  a manner  of  speaking,  the  analyses  of  Section  II 
embody  microscopic  determinism,  those  of  Section  III  embody 
macroscopic  determinism,  and  those  of  Section  IV  embody  mac- 
roscopic statistical  analysis.  It  is  indeed  reassuring  to 
see  the  agreement  achieved,  and  there  is  always  merit  in  "see- 
ing a problem  from  different  viewpoints.  However,  we  have 
so  far  considered  only  one-sided  roughness,  and  some  ques- 
tions need  to  be  answered  before  more  difficult  problems  can 
be  successfully  tackled.  For  example,  most  of  the  macroscopic 
deterministic  analyses  (refs. 11-14)  take  advantage  of  periodic- 
ity or  quasi-periodity  (high  spatial  correlat j.on)  to  derive 
load  support  relations,  whereas  the  statistical  analysis  of 
ref.  (18)  invokes  randomness  (low  correlation) . Yet  both  ap- 
proaches yield  the  same  results! 

The  answer  to  this  paradox  lies  in  the  fact  that  sta- 
tistical analyses  of  Reynolds  equation  have  depended  on  the 
two  postulates  leading  to  eqs.  (3.1)  and  (3.2),  generalized. 
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however,  to  allow  for  weak  spatial  dependence. + But  once 
these  postulates  are  accepted,  there  is  really  little  need 
for  statistics!  For  example,  generalizations  of  the  mass- 
flux  relations  (3.1)  and  (3.2)  then  follow,  whether  or  not 
the  striations  are  random.  It  is  undeniably  desirable  to  be 
aware  of  the  likely  response  of  a film  to  haphazard  surface 
fluctuations  but,  it  is  also  desirable  to  separate  the  cru- 
cial from  the  non-crucial  arguments  in  a derivation. 

Since  the  validity  of  the  Reynolds  roughness  equations 
hinges  on  the  correctness  of  the  mass-flux  postulates , one  may 
legitimately  ask  when  they  are  valid.  As  Christensen  (18) 

emphasizes,  and  as  all  investigators  agree,  in  the 

limit  of  short  roughness  wavelength.  This  conclusion  is  man- 
ifest by  the  numerical  studies  cited  in  Section  II,  and  sup- 
ported by  the  multiple-scale  analysis  in  (15)  . However,  if 
these  postulates  are  taken  as  the  starting  point  for  a rough- 
ness analysis,  all  trace  of  wavelength  disappears  from  the 
equations!  So  there  can  be  no  self-checking  for  consistency. 
Nevertheless,  these  postulates  have  physical  relevance  and 
irresistible  simplicity,  so  all  of  us  will  doubtless  continue 
to  use  them  in  new  applications,  confident  that  experiment 
and  further  mathematics  will  later  corroborate. 

Mathematically,  the  point  to  be  made  concerning  the 
postulates  is  that  they  constitute  guesses  concerning  solu- 
tions of  Reynold  equation.  If  one  is  to  use  Reynolds  equa- 
tion, why  not  use  it  from  start  to  finish?  The  multiple- 
scale  analyses  in  refs.  (15,21)  essentially  accomplish  this 
objective  for  one-sided  roughness.  The  series  in  A/L=e  of 
eq.  (3.1)  is  carried  through  to  second-order,  and  interior 


+ The  analysis  of  Tzeng  & Saibel  (16)  postulates 
that  the  total  bearing  throughput,  Q,  will  be  invariant. 
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solutions  (to  eqs.  like  (4.9)  and  (4.10))  are  asymptotically 
matched  at  film  peripheries  by  Mui j derma n- type  boundary  so- 
lutions (22) . Additional  aspects  of  film  shape  manifest 
themselves  in  the  second-order  treatment,  with  rectangular 
grooving  preserving  striking  simplicity.  The  internal  con- 
sistency of  the  whole  treatment  is  very  convincing,  but  the 
analytical  difficulty  is  such  as  to  discourage  immediate 
attempts  on  more  complicated  roughness  patterns  (two-sided, 
two-dimensional , etc . ) . 

While  difficulties  have  been  pointed  out  in  the  deri- 
vation by  statistics  of  differential  equations  for  Reynolds- 
type  roughness,  there  can  be  no  doubt  that  statistics  are 
often  necessary  to  describe  properties  of  the  roughness  it- 
self. In  order  words,  the  coefficients  in  differential  equa- 
tion such  as  (4.9)  and  (4.10)  indicate  what  properties  of 
the  local  film  contour  are  required,  and  statistical  distri- 
bution functions  such  as  (4.7)  give  our  best  estimate  of 
their  values. 

On  the  other  hand,  it  is  writer's  opinion  that,  at  the 
present  time,  much  can  be  still  be  learned  by  using  simple 
roughness  functions  and  the  short-wavelength  postulates. 
Statistics  should  be  used  sparingly,  and  principally  to  pro- 
vide numerical  estimates  of  various  roughness  parameters. 

A conceptual  difficulty  which  the  smoothening  and  statis- 
tical approaches  have  in  common  is  the  definition  of  sample 
space  for  local  averages.  Both  suppose  the  existence  of 
such  local  averages  to  be  not  critically  dependent  on  the 
choice  of  sample  size.  The  problem  here  is  precisly  the  same 
as  the  definition  of  density  for  a gas.  For  too  small  a sam- 
ple volume,  fluctuations  in  the  value  of  density  occur  as 
individual  molecules  are  included  or  excluded  from  the  volume. 
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VI.  Two-Sided  Reynolds  Roughness 

In  1971  Berthe  & Godet  (23)  pointed  out  that  roughness 
analyses  up  to  that  time  were  not  appropriate  for  films  with 
two  rough  bounding  surfaces  if  mutual  tangential  motion 
should  also  be  present.  Let  us  now  consider  situations  where 


Figure  6 . 1 Two  Rough  Moving 
Surfaces . 

such  surfaces  are  involved.  In  these  instances,  the  rough- 
ness on  the  moving  surface  "sees"  a varying  mean  film  thick- 
ness. For  example,  in  Figure  6.1,  surface  *1  is  "moving"  if 
Uj,*U2.  A steady-state  case  with  transverse  roughness  will  be 
treated.  The  f ilm-thickness  distribution  is  given  by: 

h = h(x)  + 0^(x-Uj^t)  + 02(x-U2t)  [6.1] 


and  the  applicable  Reynolds  equation  for  an  incompressible 
fluid  is: 


D h^  3p  . 3 h^  ap  _ ah  , 3h  * 

J3^TIT['3lf‘'TyT2U35“  ^ ' Tt 


* Velocities  are  restricted  to  the  x-direction 
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Substitution  of  the  film- thickness  expression  (6.1)  into 
(6.2)  yields: 


)x\ 


12n 


iB  - 
ax 


(U^+U2)^  (U^-U2) 
— = — h+  — = — 


(63^-62) 


1+  -1  ^ 
J 3y  I2p 


lE  = 

3y 


[6.3] 


For  the  case  of  two-dimensional  contact  of  cylinders 
in  pure  rolling,  U^=U2=U  and  a/ay=0,  and  eq.(6.3)  reduces  to: 


a ap  _ ^ah 
ax  i2p  ax  ^ax 


[6.4] 


This  equation  can  be  explicitly  integrated  subject  to  bound- 
ary conditions  on  "p".  Michau,  Lafont,  Berthe  and  Godet  (24) 
and  Tallian  (25)  have  used  it  it  estimate  the  pressures  and 
pressure  excursions  due  to  rippling  which  may  occur  in  roll- 
ing contact. 


When  rippling  effects  remain  in  the  convective 

term  of  Reynolds  equation.  Thus  for  the  simple,  infinitely 
wide  slider  in  Fig.  6.1,  with  U2=0,  eq.  (6.3)  becomes: 

. 0 16.51 


Chow  and  Cheng  (26)  invoke  a generalized  form  of  the  flux 
continuity  postulate,  hypothesizing 


Qx 


IE  U 
ax  2 


(h-ej^+e2) 


[6.6] 


as  locally  smooth  and  continuous.  Precisely  as  in  eq.  (3.2), 
one  then  obtains: 


Qx<73’m 
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_L-(1E) 

12y  ax  m 
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4 


[6.7] 
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With  taken  as  independent  of  "x" , Chow  & Cheng  eliminate 
it  by  differentiation  to  obtain  a second-order  differential 
equation.  This  they  apply  to  the  elastohydrodynamic  contact 
problem.  Other  forms  of  differential  equation  are  derived 
by  them  for  longitudinal  roughness,  and  for  different  combi- 
nations of  U2. 

The  Chow-Cheng  differential  equations,  although  de- 
rived in  a different  manner  (i.e.,  on  a postulatory  basis), 
are  special  cases  of  the  differential  equation  published  in 
ref.  (27) . This  earlier  derivation  was  via  multiple-scales, 
and  again  in  the  interest  of  brevity,  we  shall  forgo  an  ac- 
count of  the  method.  It  is  recommended,  however,  to  those 
who  would  like  to  contribute  new  results  in  this  field.  The 
final  result  in  (27)  is:‘*' 


_8 

3t 


(HtTq ) + ( Aj^+A2 ) • V (Hit Q ) + 


] = 


s 3 

' "o-ai' 


3n' 


[6.8] 


Here  the  coordinate  system  consists  of  "n"  normal  to  the 
striations  (^,  the  unit  normal)  and  "s"  parallel  to  the  stri- 
ations.  The  striations  on  one  surface  must  be  parallel  to 
those  on  the  opposing  surface,  but  the  fluid  may  be  compres- 
sible (iiq=1  in  all  coefficients,  if  incompressible),  "p" 
may  be  arbitrary,  and  the  flow  may  be  transient.  In  most  in- 
®1  ® 2 

stances,  I,  = (—7)  and  Io(— t)^,  as  in  eq.  (6.7). 
i m i m 


* 


"n"  and  "s"  are  here  made  non-dimensional  with  "L". 
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For  the  case  of  transverse  roughness  with  a simple 
slider  bearing,  the  effects  of  location  of  the  roughness 
are  shown  in  Fig.  6.2. 


O Z 4 G ^ fO 

3<t/c 


Figure  6.2  Load  Capacity  as 

Function  of  Transverse 
Roughness  Location 
(from  ref. 27) 

The  maximum  film  thickness  excursion  due  to  roughness  is  the 
same  for  each  curve,  whether  from  one  surface  alone,  or  two 
in  combination.  The  manner  in  which  load  capacity  varies 
with  location  can  be  rationalized.  Essentially,  pressure 
buildup  is  the  result  of  fluid  resistance  to  being  drawn 
by  convective  shear  into  a narrowing  film.  For  given  surface 
speed  and  mean  film  thickness,  less  resistance  will  be  expe- 
rienced if  some  of  the  fluid  is  drawn  along  in  the  inter- 
stices of  a moving  rough  surface.  Hence  the  lowest  load  ca- 
pacity is  associated  with  the  "rough  moving"  curve  of  Fig. 
6.2,  with  intermediate  and  highest  capacities  resulting  from 
transfer  of  part,  or  all,  of  the  roughness  to  the  stationary 
surface . 
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The  argument  advanced  above  would  lead  one  to  conjec- 
ture that  roughness  would  have  less  effect  with  a compres- 
sible fluid,  since  then  the  fluid  "has  the  choice"  of  leav- 
ing a location  by  pressure  gradient,  or  simple  being  com- 
pressed within  the  roughness  depressions.  Some  examples  in 
ref.  (28)  illustrate  the  phenomenon. 


VII.  Two-Dimensional  Reynolds  Roughness. 


The  next  logical  generalization  of  treatment  is  to 
non-striated , or  two-dimensional  Reynolds  roughness.  Natu- 
rally, as  we  advance  in  difficulty  the  applicable  literature 
thins  out. 


For  isotropic  Reynolds  roughness,  Christensen  & Tgsnder 
propose  (see  also  eq.4.11) 
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7.1 


To  this  writer's  knowledge,  no  proof  of  this  relation  has 
appeared  in  the  open  literature.  For  reasons  to  be  given 
below,  I believe  it  to  be  incorrect. 

There  appears  to  be  the  opportunity  for  further  pro- 
gress along  the  lines  of  Hamilton,  Walowit  and  Allen  (29)  and 
Tsao  and  Tong  (30).  These  investigators  determine  by  Reynolds 
equation  the  flow  about  a single  cylindrical  asperity  with 
flat  or  spherically-capped  top  (Fig.  7.1). 


Asperity  Geometry 
(from  ref.  29) 


Figure  7 . 1 


To  mimic  surface  roughness,  they  then  form  a grid  oi  such 
asperities  (Fig.  7.2). 


O O 


Figure  7.2.  Asperity  Distributions 
(from  ref.  29) 


An  externally-impressed  pressure  gradient  is  not  included  in 
their  analyses,  their  interest  being  chiefly  in  face  seals. 
Cavitation  is  included,  but  only  in  the  Gumbel  sense  with 
"lopping  off"  of  negative  pressures. 

For  a single  asperity  of  the  form  shown  in  Fig.  6.1, 
where  the  film  thickness  has  two  different  constant  values, 
Reynolds  equation  reduces  to  Laplace's  equation 

v2p  = 0 [7.2] 

in  both  inner  (r<RQ)  and  outer  (r>RQ)  regions.  Suitable 
solutions  are; 

Inner:  p = r 

} I V . n 

( * t • 

t )i  1 1 < • I : I ' k I I I k , / I 


Matching  mass  fluxes  and  pressures  at 
interface  (r=RQ) , we  obtain; 


the  regional 
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Inner: 

P “ 

k®x 

KX 

Outer: 

P = 

k®x 

+ 

Kx(^)2 

where : 

K » 

(a+b) 

^-b^}-6wUa 

(a+b)Hb^ 


[7.5J 


The  above  pressure  distributions  consist  of  an  sxternally- 
Inpressed  gradient  and  a deviational  component.  Typical 
results  for  the  latter  are  shown  In  Fig.  7.3. 


Figure  7.3.  Typical  Deviational  Pressure 
Distribution  (O^O)  and 
Isobars  for  Single  Asperity 
(from  ref. 29) 


Now  although  the  disturbance  (deviational)  pressure  vanishes 
at  r -*  •/  there  remains  a net  ‘•drag"  resulting  from  the  as- 
perity. Thus,  limit  M 

Drag  - [ I pdyl  - [ f pdyl  - j(kj2L)dy  17.6] 
[-M  Jx— L L-M 
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The  result,  independent  of  L>Rq,  is; 

[12vUa-2k®{  (a+b)  ^-bM] 

Drag  = [7.7] 

(a+b) 3+b3 

No  net  disturbance  effect  exists  for  the  flow. 

Now  if  many  protuberances  were  present  on  a surface, 
as  a first  approximation  one  might  consider  each  immersed 
in  an  external  pressure  gradient,  partly  impressed  by  outside 
agencies  and  partly  an  integrated  effect  from  other  asperi- 
ties. A macroscopic  theory  relating  average  flow  to  surface 
velocity  and  average  pressure  gradient  might  thereby  be 
built.  This  analysis  has  yet  to  be  done. 

For  the  case  of  purely  Poiseuille  flow  (Uj^=U2=0)  a 
satisfactory  analysis  can  be  made.  The  differential  equa- 
tion is: 
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3x 
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-3  h3|£  = 0 
3y  8y 


(7'. 8] 


This  is  the  heat  conduction  equation  with  "p"  playing  the  role 
of  temperature, "h3 " playing  the  rSle  of  thermal  conductivity 
and  lineal  mass  flow  playing  the  rSle  of  heat  flux. 


In  ref. (31)  it  is  proved  that  the  heat  flow  through  a 
region  will  be  less  than  that  with  arbitrarily  assumed  shapes 
for  the  isotherms,  and  greater  than  that  with  arbitrarily 
assumed  lines  of  heat  flow  (adiabats) . Translated  to  Rey- 
nolds equation,  this  statement  becomes:  the  lineal  mass  flow 
will  lie  less  than  with  arbitrarily  assumed  Isobars  and  gr«^uter 
than  that  with  arbitrarily  assumed  streamlines.  Let  us  apply 
this  theorem  to  the  roughness  pattern  shown  in  Figure  (7.4). 
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Figure  7.4  Two-Dimensional  Roughness 
Pattern 


First,  assume  that  the  isobars  are  horizontal.  Then  for  the 
strip  "6y"  lying  along  A-A,  "fip"  must  be  constant.  The  local 
mass  flow  is 


The  integrated  mass  flow  along  the  stretch  "W  is: 

W 

(6y)My  = - Y5^|h3dxdp  = - 


[7.9] 


Or: 


12  vM  ^ 

dp  = - — 


[7.10] 


[7.11] 


A second  integration,  this  time  with  respect  to  "y", 


yields : 
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p(L)-p(0)  = - 
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[7.12] 


Now  so  long  as  "W"  and  "L"  contain  many  roughness  elements, 
eq.  (7.12)  is,  effectively: 
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and  this  is  an  upper  limit  for  the  flow. 

Alternatively,  one  can  assume  vertical  streamlines, 
such  as  those  for  "6x"  between  B-B  in  Figure  '(7.4).  In  this 
case,  for  each  element  "6x": 


(m  )h"^  = - 
y 12y  3y 

• t -3  * 

p(L)-p(0)  = - 12p(my)!h  ■’dy  = - 12ymyh  ^ L 


[7.14] 

[7.15] 


With  the  gross  pressure  difference  p(L)-p{0)  treated  as  con- 
stant, integration  is  now  performed  along  "W"  to  give: 


M 


_ _ 1 p(L) -p(0) / 1 


12p 


[7.16] 


A lower  limit  for  the  flow  is  thus  obtained: 

(m,.)„  = - 


y m 


""  ii 


3y 


[7.]7] 


f 
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Potential  error  in  the  estimate  for  (m  ) is  minimized 

y m 

by  use  of  the  arithmetic  average  from  eqs.  t7.13)  and  (7.17) 
Thus : 


[7.18] 


For  the  case  of  heat  conduction  and  a configuration 
similar  to  that  shown  in  Fig. (7.4),  Yovanovich (32)  has  found 
the  error  of  the  equivalent  of  eq.  (7.18)  not  to  exceed  5%. 

For  isotropic  roughness  the  indicated  averaging  processes  must 
be  independent  of  direction.  Equation  (7.18)  then  reduces 
to: 


+ h^lvp 


[7.19] 


where  the  indicated  averages  are  over  areas  local  to  the  point 
under  consideration.  Equation  (7.19)  differs  from  (7.1)  and 
this  writer  believes  eq.(7.19)  to  be  more  accurate. 

Equation  (7.18)  is  actually  valid  for  two-sided  rough- 
ness, but  has  only  been  deduced  for  the  case  of  no  surface 
motion.  Obviously,  the  subject  of  two-dimensional  roughness 
requires  further  research. 


All  analyses  so  far  discussed  have  been  for  "Reynolds 
roughness"  and  have  involved  the  assumption  of  the  adequacy 
of  Reynolds  equation  to  describe  the  flow  within  the  film. 

Now  two  conditions  have  to  be  met  in  order  that  Reynolds 
equation  be  a proper  approximation  to  the  Navier-Stokes  equa- 
tions. First,  the  film  Reynolas  number,  pU^h/v,  must  be  small. 
Second,  the  change  in  film  thickness  per  unit  of  film  thick- 
ness must  be  small  (i.e.,  dh/dx<<l) . Whereas  the  first  of 
these  conditions  is  very  often  satisfied,  the  second  condi- 
tion may  be  violated  by  roughness  contours,  even  when  the 
slope  of  mean  film-thickness  is  very  smal} . In  such  instances, 
Stokes  equations,  not  Reynolds,  must  be  used. 


In  1962,  Citron  (33)  used  Stokes  equations  to  calcu- 
late the  effects  of  longitudinal  roughness  in  the  gap  between 
concentric  rotating  cylinders.  But  the  effects  of  transverse 
roughness  are  more  interesting.  In  an  recent  paper.  Sun  & 
Chen  (34)  have  investigated  transverse  Stokes  roughness  ef- 
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and  the  slider  bearing  in  Fig.  4.1  will  serve  for  nomencla- 
ture . 


In  the  limit  of  small  "e",  the  terms  on  the  rhs ' s of 
eqs.  (8.2)  and  (8.3)  disappear,  and  Reynolds  relations  are 
recovered.  Sun  & Chen  solve  these  simpler  relations  subject 
to  actual  boundary  conditions,  including  roughness,  and  insert 
the  results  into  the  rhs's  to  obtain  improved  solutions 
("first  iterations").  Figure  8.1  shows  their  results  for 
load  capacity  of  a flat  slider  bearing. 


Figure  8.1  Load  Capacity  as  Function 
of  Random  Roughness  for 
Simple  Slider  (h  /hmin=2*7) 

IT\^X 

Stolces  & Reynolds  Roughness 
Theories  Compared  (from  ref. 
34)  . 

The  presentation  in  Fig.  8.1  is  intended  to  reflect 
variations  about  an  operating  condition  of  L'\-0(25mm),  c'v-O 
(12.5nm),  av0(*25vim),  and  A'\-0  ( 12  . Spm)  . Under  these  circum- 
stances, Reynolds  roughness  theory  predicts  virtually  no 
effect,  whereas  the  authors'  theory  predicts  as  much  as  30% 
load  increase  at  0/c=O.O6  and  a roughness  wavelength  equal 
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to  the  minimum  film  thickness  (i.e.,  c/A=l) I On  the  other 
hand,  when  c/a=1/2tt,  the  predicted  load  increment  only  slight- 
ly exceeds  that  obtained  by  Reynolds  equation.  The  curves 
are  for  typical  random  roughness. 

In  his  doctoral  dissertation,  Rhow  (28)  studied  the 
effects  of  transverse  and  longitudinal  rectangular  grooving 
by  means  of  accurate  numerical  solutions  of  the  biharmonic 
Stokes  equation. 


y ' 

f K 


Figure  8.2  Nomenclature  for  Rectangular 
Grooving . 


This  particular  roughness  configuration  (see  Fig.  8.2)  was 
chosen  because  of  its  simplicity  for  computer  calculations, 
its  relevance  to  bearings  with  purposeful  grooving,  and  its 
indicative  possibilities  for  other  types  of  roughness.  The 
lineal  flowrate  corresponding  to  (a)  U*0 , Ap=0  and  (b)  U=0, 
Ap*0  over  cne  wavelength  was  calculated,  and  presented  in  the 
following  form  (for  transverse  roughness) ; 


m 


b\^  = ^ 
[l2\i  ^-3  3x. 


[8.5] 


The  coefficients  "A"  and  "B"  are  unity  for  h/A  -»  oo( Reynolds 
roughness),  and  their  deviation  from  unity  is  thcrolere  a 


39 


measure  of  the  importance  of  using  Stokes  equation. 

Over  wide  ranges  of  geometric  parameters  (d/a, 
a/A)  the  coefficient  "A"  was  found  to  lie  fairly  near  unity, 
indicating  that  Reynolds  roughness  results  for  such  flow 
will  be  reasonably  accurate.  On  the  other  hand,  as  shown  in 
Fig.  (8.3)  the  Poiseuille  coefficient  "B"  can  drop  substanti- 
ally. 


Figure  8,3  Poiseuille  Coeff,"B"  for 

Rectangular  Grooving  (a=A/2) 
(adapted  from  ref.  28) . 

To  use  the  above  results  in  a bearing  calculation,  one 
invokes  the  short-wavelength  flux  postulate,  and  uses  local 
values  of  the  film  constructs  A,  B,  h”'/ , h"3 , etc.,  in  eq. 
(8.5).  This  equation  is  then  treated  as  a differential  equa- 
tion with  liij^  constant.  Figure  (8.4)  shows  Rhow's  results* 
for  a flat  slider.  The  difference  between  Stokes  & Reynolds 

+ For  these  calculations , "A"  was  taken  as  unity,  and 
"B"  was  evaluated  in  a manner  consistent  with  Fig.  8.3. 
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results  is  here  much  less  than  in  the  work  of  Sun  & Chen 
(Fig. 8 . 1)  . 


o -z  ■4_  ■&  S JO 


Figure  8.4  Load  Capacity  as  Function  of 
Roughness  (rectangular  groov- 
ing)  for  Simple  Slider 
(h„,^/h  . =1)  Stokes  Analysis 
(from  ref.  28) . 

When  reasonable  allowances  are  made  for  differences 
in  roughness  configuration,  bearing  film  convergence,  etc., 
the  differences  in  the  results  between  Figs.  8.1  and  8.4 
remain  inexplicably  large.  Furthermore,  the  Sun-Chen  con- 
clusions tend  to  restrict  drastically  the  applicability  of 
Reynolds  roughness  theory.  Accordingly,  this  writer  has 
undertaken  a different  analysis  which  is  less  vulnerable  to 
criticism  than  those  previously  made.  The  gist  of  this  analy- 
sis, and  its  conclusion  will  here  be  presented. 

Invoke  the  short-wavelength  hypothesis,  and  consider 
two  parallel  plates,  the  upper  one  moving,  the  lower  one 
rough,  as  in  Figure  8.5.  In  Stokes  approximation,  the  stream 
function  satisfies  the  biharmonic  equation. 
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Thus : 


V**!!;  = 0 


[8.6] 


Both  externally  impressed  shear  and  pressure  gradient  are 
accounted  for  by  the  following  boundary  conditions: 


i|((x,ef(x))  = (x,ef  (x)  ) = 0 
ij)y(x,h)  = U ; \|»(x,h)  = ^ 

An  expansion  is  made  in  terms  of  amplitude,  "6", 
so  that: 


[8.7] 


i|/  = i]<Q(x,y)+6i(ij^  (x,y)+62ij^2  +• 


[8.8] 


With  this  expansion,  there  is  always  a "6"  small  enough  to 
accommodate  any  chosen  roughness  wavelength, A . The  analysis 
must  be  carried  through  to  second  order,  since  first-order 
effects  on  load  vanish. 

Tt  it  both  impossible  and  1 nnppi  opr  laU*  lo  (jivo  de- 
tails of  the  analysis  here, but  tlic  conclusions  can  be  summa- 
rized as  follows: 
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In  eq.  (8.5):  A = l-6^f^(F*4) 


B = l-362f7(F-4) 


where  6 = (roughness  amplitude) /h 


f(x)  = roughness  shape  function 


expressed  as  = ZB^sinw^x  in  Fourier  Series. 


f(x)Z  F(hw  . )Bj^sin(a)^x) 

F*  = — — 

f^TxT 

and:  F(C)  = [tanh£{ l+gtanhg }-C ] 

(1+e^) tanh^s-e^ 


[8.9] 


[8.10] 


[8.11] 


[8.12] 


Equations  (8.9-8.12)  were  used  in  eq.  (8.5)  for  the 
slider  bearing  of  Fig.  8.1,  with  = 0.06,  and  with 

sinusoidal  roughness. 


Orel  34)  ^ 


^ ^ 

ds 


Z 4-^,6  i to 

T,M 

Figure  8.6  Load  Capacity  as  Function 

of  Sinusoidal  Roughness  for 

Simple  Slider  (^  . _=2  • 7) 

,c  A max  min 

o/h  . -0.06. 

' min 
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The  resulting  curves  are  shown  in  Fig.  8.6.  Again  the  pre- 
dicted magnitude  of  the  Stokes  roughness  effects  differs 
markedly  from  the  corresponding  prediction  of  Sun  & Chen. 
Since  I am  familiar  with  the  ingredients  of  two  independent 
theories  which  conflict  with  the  Sun-Chen  results,  I have 
to  conclude  the  Sun-Chen  results  are  wrong.  Why,  I am  not 
certain,  but  possibly  the  iteration  process  they  adopt  does 
not  converge  for  A'vh.  For  example,  in  general: 


so  that  in  eq.  (8.2) 


3^u' 


2 


Thus  the  rhs  of  (8.2)  is  not  small  compared  with  the  Ihs 
when  h/A=l.  No  doubt,  with  this  questioning  of  their  results, 
the  authors  will  recheck  their  work  and  we  shall  hear  more 
of  the  matter. 
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There  are  occasions  when  it  is  pointless  to  consider 
roughness  effects  without  simultaneously  considering  mole- 
cular slip  effects.  For  example,  the  "flying  heads"  used 
for  information  retrieval  from  rotating  disks  in  modern 
computers  are  actually  small  air  slider  bearings.  Minimum 
mean  clearance  in  these  devices  is  now  being  reduced  to  ap- 
proximately 0.4pm.  At  the  same  time,  the  mean-free  path  of 
the  air  molecules  is  about  O.lym  and  the  roughness  amplitude 
of  the  surface  finish  is  about  O.Oipm. 


The  first  effects  of  rarefaction  are  customarily  ac- 
counted for  by  the  use  of  slip-flow  aerodynamics.  Thus, 
within  the  flow  interior,  standard  continuum-fluid  differen- 
tial equations  are  employed,  while  failure  of  the  macroscopic 
fluid  velocity  to  agree  with  that  of  its  bounding  surfaces 
is  modelled  by  the  Maxwell  slip  condition.  Thus: 

Au  = a|^  [9.1] 


Here  "Au"  is  the  above-mentioned  tangential  velocity  dis- 
crepancy, "a"  is  the  mean-free  molecular  path,  and  "3u/3n" 
is  the  velocity  derivative  normal  to  the  surface.  When  eq. 
(9.1)  is  used  in  connection  with  Reynolds  equation  (Burgdor- 
fer,  1959  ,’^)  the  lineal  mass  flow  is  given  by  the  following 
simple  equation: 


Uh 

P 2 


[9.2] 


In  the  presence  of  ono-sJded  roughness,  smoothoning 
by  means  of  the  flux  continuity  postulate  and  the  load  in- 
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coinpressibity  assumption’"  yields,  to'  r^i  /ecie  roughness: 


where : 


12yRTLm  ^-2  , 

X _ „ . 1 i.T 

2 3 "^-3  ’^o  ' 

p^c-^  X ^ J O 

3. 


[H^  + 6—  — ] ^ 
^ ^o 


[9.3] 


[H^+6—  ~]  ^ 

C IT 

o 


Correspondingly,  for  longitudinal  rougl  ess,  one  gets; 


12yRTL*  , S „3  3tto 

~ — ^in  = A Hu  -H  IT  

2^3  y y o o3n 

^a 


"^3  3 T4^ 

where:  H = (H  +6- — — ) 

C IT 

o 


[9.4] 


[9.5] 


[9.6] 


The  components  can  be  combined,  as 
lineal  mass  flux  vector. 


•=q  ' 3,  , to  form  the 


Now  the  foregoing  analysis  brea  .s  down  for  short 
roughness  wavelengths.  Recently  Chang  and  this  author 
have  performed  kinetic-theory  calculati ■ .is  for  rectangular 
grooving.  In  this  work  the  Boltzmann  equation  is  solved  by 
dividing  velocity  space  into  iour  i.so-energetic  groups  of 
different  direction.  The  method  is  suitable  for  low  macro- 
scopic speeds,  and,  in  appropriate  l-imits,  "Ipids  accurate 
slip  flow  or  continuum  results,  .on;_  .;ete  overlap 

with  the  continuum  work  of  Rhow  and  Elrod  i;:  obtained.  It  is 


* Both  are  endorsed  by  an  unpub  .i. shed  multiple-scale 


treatment  of  the  author. 
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hoped  that  details  will  be  published,  but  the  gist  of  the 
principal  conclusion  can  be  stated  as  follow: 

"When  the  effects  of  molecular  slip  and  of  one-sided,  trans- 
verse Stokes  roughness  are  simultaneously  present,  one  can  use 
eg.  (9.3),  modified  by  the  same  continuum  correction  factors 
as  would  be  used  with  eq.(8.5)  for  the  same  roughness  contour" 
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X.  Conclusions  and  Acknowledgment 

Although  the  accomplishments  of  analysis  in  predicting 
the  effects  of  roughness  on  laminar  lubricating  films  are 
perhaps  encouraging, the  work  is  incomplete,  and  a number  of 
differences  in  outlook  and  results  remain  to  be  reconciled. 

Two-sided  Reynolds  roughness  has  been  analyzed  only 
for  parallel  striations  on  opposing  surfaces.  Two-dimension- 
al (i.e.,  non-striated) , roughness  has  been  treated  only 
for  stationary  surfaces.  Even  the  first  effects  of  Stokes 

roughness  are  in  some  dispute,  though  there  is  agreement  that  | 

{ 

Reynolds  equation  is  not  as  suitable  as  we  might  like  for  | 

roughness  analyses, due  to  its  long-wavelength  restriction. 

Much  of  the  more  recent  work  would  benefit  from  a less  heu- 
ristic development.  And  finally,  the  areas  for  utility  of 
statistical  analyses  need  to  be  better  delineated. 

The  writer  is  most  grateful  to  his  hosts  at  the  Tech- 
nical University  of  Denmark  for  providing  time  and  facilities 
for  the  preparation  of  this  paper.  A basis  for  its  contents 
was  provided  by  experience  and  results  acquired  during  re- 
search contracted  by  the  U.S.  Office  of  Naval  Research  under 
Task  Nr.  062-360,  N00014-67-A-0108-0021 , and  conducted  at 


Columbia  University,  New  York,  10027. 
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NOMENCLATURE 

(Some  symbols  which  are  used  in  a restricted  portion  of  the 
text  are  defined  there  on-the-spot) . 

A correction  coefficient  in  Reynolds  eq. , eq.8.5 
B correction  coefficient  in  Reynolds  eq.,  eq.8.5 
c nominal  clearance  of  bearing  (=  min.  cl.  for  sliders) 
parameter  in  roughness  distribution,  eq.4.7 
C -C^,  friction  force/load 

h film  thickness,  including  roughness  h,  mean  local  film 
thickness,  h = h+e  ; h^  = specific  value  of  0 
H h/c 

L characteristic  gross  dimension  of  bearing 
m slope  of  slider  bearing  t^jj,ax”*^min^ 

m,  lineal  mass  flow  rate  per  unit  transverse  distance 
n coordinate  perpendicular  to  striations 

Also,  number  of  waves  used  in  deterministic  calculations 
p fluid  pressure 

Q lineal  volumetric  flow  in  film.  Volume  per  unit  time  per 
unit  distance  transverse  to  flow 
r distance  from  asperity  center 
R R^  asperity  radius 
s coordinate  parallel  to  striations 
t time 

u fluid  velocity,  x-direction 
U surface  velocity  in  x-direction 
v fluid  velocity,  y-direction 
W (load)c^/uUL^ 

X cartesian  coordinate 
y cartesian  coordinate 

z coordinate  perpendicular  to  striations.  See  Fig.  1.1 
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e orientation  angle  of  striations.  See  Fig.  1.1 
6 roughness  amplitude/mean  film  thickness 
A striation  wavelength 
e A/L  or  c/L 
t z/A 

n y/L  also,  fluid  viscosity 

e circumferential  angle,  or  magnitude  of  roughness 
X mean-free  path  of  molecules;  X = ambient  X 
A 6viUL/P  c ; for  liquids,  use  p = P(-3^)m 

Si  3 V p 1 

P fluid  viscosity 

o standard  deviation  of  roughness 
T p c^t/12pL^;  for  liquids,  use  p =P(-^)„ 

<i>  stream  function.  9iJ;/8y  = u;  3'('/9x  = -v 


Special  Symbols  & Subscripts  ; 

m mean  or  local  space  average  . 

7 T mean,  statistical  or  space  average 

< > expected  value,  or  statistical  mean  ^ 

subscript  zeroth  order,  first  order,  etc.  function  of  an 
0,1,2,  expansion,  or  used  to  signify  location, 
a ambient 

- i 

n unit  normal  in  n-direction  j 

- ■ i 

j unit  normal  in  s-direction  t 


j 

1 
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